function [ma] = SolveForMA(L, tau, theta, Counties)

    x0 = ones(1,Counties);
    Tol = 10^-5;
    new = FunctionForMA(x0,tau,L,theta,Counties);
    x = x0;

    while (sqrt(sum((new - x).^2)) > Tol)
        x = new;
        new = FunctionForMA(new,tau,L,theta,Counties);
    end

    ma = new.^-theta;

end
